Chandra Observations of Diffuse X-Rays from the Sagittarius B2 

Cloud 



O 

o 

(N 



Hiroshi Murakami, Katsuji Koyama 
Department of Physics, Faculty of Science, Kyoto University, Sakyo-ku, Kyoto 606-8502, Japan; 
hiro@cr.scphys.kyoto-u.ac.jp, koyama@cr.scphys.kyoto-u.ac.jp 

and 

Yoshitomo Maeda^ 

Department of Astronomy and Astrophysics, The Pennsylvania State University, 525 Davey Lab. 
University park, PA 16802; maeda@astro.psu.edu 



> 

m 
(N 

in 
o 



o 



ABSTRACT 

We present the first Chandra results of the giant molecular cloud Sagittarius B2 (Sgr B2), 
located about 100 pc away from the Galactic center. Diffuse X-rays are clearly separated from 
one-and-a-half dozen resolved point sources. The X-ray spectrum exhibits pronounced iron K- 
shell transition lines at 6.40 keV (Ka) and 7.06 keV (K/?), deep iron K-edge at 7.11 keV and 
large photo-electric absorption at low energy. The absorption-corrected X-ray luminosity is ~ 
1 X 10^^ erg s~^, two orders of magnitude larger than the integrated luminosity of all the resolved 
point sources. The diffuse X-rays come mainly from the south-west half of the cloud with a 
concave-shape pointing to the Galactic center direction. These results strongly support the 
ASCA model that Sgr B2 is irradiated by an X-ray source at the Galactic center side. 



Subject headings: Galaxies: Milky Way 
Individuals Sgr B2 

Introduction 



Interstellar: clouds — X-rays: sources — X-rays: spectra 



The central region of our Galaxy (here, GC) 
is behind the large extinction in the optical band, 
but visible with the radio and infrared bands. The 
GC has been revealed to exhibit highly complex 
features, which are possibly attributable to either 
a putative massive black hole (MBH) of '--^ 3 x 10^ 
Mq at the position of the radio non-thermal source 
Sgr A* (Ghez et al. 2000) or star burst activity in 
the GC. 

The interstellar gas to the GC also becomes 
transparent at higher energy X-rays. Diffuse X- 
ray emission in the GC was discovered with the 
Einstein Observatory (Watson et al. 1981). The 
Ginga satellite then found Ka lines of highly ion- 
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ized irons associated with a high temperature 
plasma in the central ^100 pc region, the ex- 
plosion energy of which is as large as lO^'* ergs 
(Koyama et al. 1989; Yamauchi et al. 1990). The 
ASCA satellite, with the imaging capability in 
the wide energy X-ray band (0.5-10 keV) cou- 
pled to the reasonable energy resolution, has fur- 
ther found a hot plasma inside the Sgr A ra- 
dio shell with an oval-shape of ~ 2' x 3' across 
(Koyama et al. 1996). Early results with the 
Chandra X-ray Observatory resolved out numer- 
ous X-ray structures from the central 20 pc (8') 
region of the GC, and found a weak point source 
of Lx — 10^"^ ergs s^^ at the position of Sgr A* 
(Baganoff et al. 2001). This is nearly one billion 
times lower than the Eddington luminosity for a 
MBH of - 3 X 106Mq. 

This "X-ray quiet" Sgr A* is in a sharp con- 
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trast to the high X-ray activity of the surround- 
ing diffuse hot plasmas. A hint to connect these 
X-ray phenomena, the quiet Sgr A* and active en- 
vironment, was obtained with the ASCA obser- 
vations on Sgr B2, a giant molecular cloud lo- 
cated ~100 pc away from Sgr A*. Koyama et 
al. (1996) discovered diffuse emission of a Ka 

6.4 keV line of neutral or low ionized irons from 
Sgr B2. The X-ray spectrum and morphology of 
Sgr B2 arc well explained by an "X-ray Reflec- 
tion Nebula (XRN)" model; the X-rays are due 
to reflection, photo-electric absorption and fluo- 
rescence from iron atoms (Murakami et al. 2000). 
However, no adequately bright X-ray source was 
found in the vicinity of Sgr B2 to fully account 
for the diffuse X-ray flux. One likely scenario, 
despite its considerable distance from Sgr B2, is 
a past activity of Sgr A*; an X-ray outburst of 
Lx ^ 3 X lO''^ ergs possibly caused by a 
surge of accretion onto the MBH in the near 
past, and is currently in a quiescent accretion 
(Koyama et al. 1996; Murakami et al. 2000). 

The limited spatial resolution of ASCA, how- 
ever, could not exclude possible contamination of 
many point sources such as young stellar objects 
in Sgr B2, which might deform the spectrum and 
morphology of the diffuse emission. We hence 
performed a high resolution imaging spectroscopy 
on the Sgr B2 region with Chandra. We suc- 
cessfully detected diffuse X-rays, resolving many 
point sources in Sgr B2. In this paper, we fo- 
cus on the diffuse X-rays, while the results of the 
point sources will be given in a separate paper. 
Throughout this paper, the distance to Sgr B2 is 
assumed to be 8.5 kpc, the same to Sgr A*, and 
is within the error of estimated distance of 7.1 ± 

1.5 kpc (Reid et al. 1988). 

2. Observation and Data Reduction 

The Chandra deep exposure observation on 
Sgr B2 was carried out on 29-30 March 2000. 
The satellite and instrument are described by 
Weisskopf et al. (1996) and Garmirc et al. (2001) 
respectively. Sgr B2 lies near the cicmtcr of the 
ACIS imaging array (ACIS-I) of four front-side 
illuminated CCDs, each with 1024 x 1024 array 
of 0.5" X 0.5" pixels covering 8'. 4 x 8'. 4 field of 
view. Data acquisition from ACIS was made in 
Timed-Exposure Faint mode with chip readouts 



of every 3.24 s. 

The ACIS front-side illuminated CCDs have 
suffered from an increase of the parallel Charge 
Transfer Inefficiency (CTI), which decreases the 
detection efficiency, energy gain and resolution. 
Sgr B2 is located at the top of amplifier 2-3 on 
chip 10, a location to be heavily affected by the 
CTI. We therefore correct these effects using the 
software described by Townsley et al. (2000), then 
no significant variation is found in the energy of 
the instrumental Ni-Ka (7.5 keV) and Au-La (9.7 
keV) lines (fluorescence due to particle bombard- 
ment on the satellite metals) at the chip loca- 
tion of Sgr B2 within < 0.5% (90% confidence). 
The energy resolution, however still varies from 
row to row: 160 and 250 eV (FWHM) at the 
bottom and top of the chip 10, respectively. In 
the spectral analysis, we use a response matrix 
based on a nearly contemporaneous observation 
of reference lines from an on-board calibration 
source (OBSID=62097), which are analyzed with 
the same CTI corrector. The effective area of the 
telescope mirrors and the detection efficiency of 
ACIS are calculated with the mkarf program in 
the Chandra Interactive Analysis of Observations 
Software (CIAO, Version 1.0). To remove back- 
ground events, we selected the ASCA grades 0, 

2, 3, 4, & 6. We also removed events from flar- 
ing pixels and artificial stripes caused probably by 
hot pixels in the frame-store region and by parti- 
cles which hit on the CCD node boundaries. 

The total count rate over the chip is increased 
by 30% in the middle of the observation due to 
a high background rate caused by the impact of 
solar energetic particles on the CCD chips. How- 
ever, we do not exclude the data during the high 
background rate because it does not damage our 
imaging and spectroscopic studies. The effective 
exposure time is 100 ks. 

3. Results and Analysis 
3.1. Point Sources 

We first search foreground X-ray sources using 
the soft (0.5-1.5 keV) band image. In the ACIS- 
I field of 17'.4 X 17'.4, we find four bright X-ray 
sources, which have Tycho-2 counterparts. The 
positions of the Tycho-2 sources have accuracy 
of 60 mas (H0g et al. 2000), hence the Chandra 
frame is fine-tuned using these four sources within 
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an absolute error of 0.25 arcsec (1 sigma disper- 
sion). 

To highlight X-ray sources near or beyond the 
Sgr B2 distance, we use the 3.0-8.0 kcV band im- 
age. The X-ray image near at Sgr B2 is compli- 
cated with diffuse emission and many faint point 
sources. For the point-source search, we exe- 
cute the CIAO 'wavdetect' software of a wavelet 
method (Freeman et al. 2000). The threshold sig- 
nificances of the wavdetect are set at 10~^ and 
0.001 for the source list and the background es- 
timation, respectively. The wavelet scales used 
are 1, v^, 2, 2v^, 4, 4v^, 8, 8v^, and 16 
pixels. We then resolve about one-and-a-half 
dozen point sources in the 3.0-8.0 keV band from 
the central 3' x 3'.5 region of the Sgr B2 cloud. 
Neither radio nor infrared counterpart is found 
from these sources, except possible association of 
the two brightest sources to a compact HII com- 
plex near the center of the cloud, Sgr B2 Main 
(de Vicente et al. 2000). 

X-ray photons are extracted from a circle of 
a radius of six times of the half power diameter 
(HPD) of point spread function (PSF) for each 
source. Since the photon numbers of individual 
sources are limited, we collect all the events from 
the point sources and make the combined point- 
source spectrum, where the background spectrum 
is taken from the point-source-free region in the 
Sgr B2 cloud shown in Figure 1. The background- 
subtracted combined spectrum is well fitted with 
an optically thin thermal plasma model with the 
best- fit parameters shown in Table 1. The inte- 
grated flux of all the point sources is 1.5 x 10~^ ph 
s"""^ cm~^ in the 2-10 keV band, which is about 8% 
of the total diffuse flux (see §3.2). The mean A^h 
value is estimated to be 1.9x10^'^ H cm^^, larger 
than that of the Galactic absorption to Sgr B2 
of ~ 1 X 10^3 H cm-2 (Sakano et al. 2000), but 
smaller than that through the Sgr B2 cloud of > 
8x10^"^ H cm^^ in the central 2' radius region (Lis 
& Goldsmith 1989, see also section 3.2). There- 
fore, we assume that most of the point sources are 
located near or in the cloud (distance is 8.5 kpc) 
and estimate the integrated luminosity (absorp- 
tion corrected) in the 2-10 keV band to be 3x 10^^ 
erg s~^. Further details of individual point sources 
will be given in a separate paper. 



3.2. DiflFuse X-rays 

Since ASCA already found the strong 6.4 keV 
line from Sgr B2, we make the 6.4 keV-line (6.15- 
6.55 keV) image to highlight the most detailed 
structure. Using the procedure given in Baganoff 
et al. (2001) and the algorithm of Ebeling, White 
& Rangarajan (2000), the image is flat-fielded, and 
is smoothed adaptively with circular Gaussian ker- 
nels. Figure 1 is the 6.4 keV line image overlaid 
on the radio intensity contours of the molecule 
"CO (Sato et al. 2000). The global morphology 
of the 6.4-keV emission is a concave-shape located 
at the south-west side (right hand side) of the 
molecular cloud pointing to the Galactic center 
direction (right hand side). The positions of the 
resolved point sources mentioned in the previous 
section are indicated by solid circles with radii of 
six times of the HPD of PSF. No significant ex- 
cess of the point sources is found in this narrow 
band (6.15-6.55 keV) image, which indicates that 
the 6.4 keV line of normal point sources, unlike 
the diffuse emission, is not prominent. The lo- 
cal excess near the center and upper-left of the 
molecular cloud are possibly due to HII complexes, 
Sgr B2 Main at (/, b) ~ (0.67, -0.04) and North at 
{l,b) ~ (0.68,-0.03). These two local excess oc- 
cupy only 5 % flux in the total diffuse flux. 

The X-ray spectrum of the diffuse emission is 
made using the data from the solid square and sub- 
tracting background taken from the dotted circle 
in Figure 1. We then fit the spectrum with a phe- 
nomenological model, a power-low function of pho- 
ton index 2.0, adding Ka and K/3 lines of neutral 
iron and absorption with the Balucihska-Church 
& McCammon (1992) cross-sections for the so- 
lar abundances (Feldman 1992). Iron abundance, 
which is sensitive to the edge depth, is treated sep- 
arately as a free parameter. The center energy and 
flux of iron Ka line are free parameters, but those 
of K/3 are linked to the Ka parameters by fixing 
the ratio to the laboratory values. Point source 
contamination is taken into account by adding the 
best-fit model of the combined point-source spec- 
trum (see §3.1). The best-fit parameters are given 
in Table 2. We find that the iron Ka line energy 
is 6.38lo keV (errors adopted here and after are 
90% confidence level unless otherwise noted), con- 
sistent with neutral irons of 6.40 keV. In order to 
examine whether highly ionized irons are present 
or not, we add 6.70 keV line, Ka of He-like iron, 
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and constrain that the equivalent width is less 
than 0.07 keV, or at most 3 % of that of neutral 
irons. Column density of neutral iron (iVpe) deter- 
mined from the K-edge structure is ^ 3.4 x 10^^ 
cm~^, while A^h determined from the low energy 
cutoff is ~ 8.8 X 10^^ cm~^. Both are nearly ten 
times larger than those of the Galactic interstel- 
lar absorption, hence are attributable to the gas 
in Sgr B2. The ASCA results, in spite of its poor 
spatial resolution, arc roughly consistent with the 
Chandra best-fit parameters, supporting that the 
contribution of the point sources is almost negli- 
gible (Murakami et al. 2000). 

4. Origin of the Diffuse X-rays 

Koyama et al. (1996) and Murakami et al. 
(2000) have argued that the diffuse X-rays of 
Sgr B2 can be explained by an XRN model; the 
Sgr B2 cloud is irradiated by strong X-rays from 
an external source, absorbs the X-rays, then re- 
flects (Thomson scattering) and/or re-emits fluo- 
rescent Fe lines. The discovery of diffuse 6.4 keV 
line from other giant molecular clouds near the Ra- 
dio Arc (Koyama et al. 1996) and at Sgr C (Mu- 
rakami et al. 2001) may indicate that the presence 
of XRNs at GC is common feature. 

The present Chandra observation clearly detect 
diffuse X-rays from Sgr B2 by resolving underly- 
ing point sources. The diffuse X-rays are emitting 
from the south-west half (right side in Figure 1) of 
the cloud with a concave-shape morphology point- 
ing to the Galactic center direction. The spectrum 
is characterized with the pronounced iron Ka line 
at 6.38lg;Q^ keV, heavy absorption at iron K-edge 
and at low energy, which confirm the ASCA results 
(Koyama et al. 1996; Murakami et al. 2000). Fur- 
thermore, K/3 line at 7.04lg;Q^ keV is also discov- 
ered for the first time. These results give further 
confirmation of the XRN concept of Sgr B2. We, 
however, examine more critically for the origin of 
the diffuse X-rays based on the new observational 
facts. 

4.1. Young Stellar Objects 

The giant molecular cloud Sgr B2 is known 
to be one of the richest regions of ongoing star 
formation in the Galaxy (e.g., de Vicente et al. 
2000). In fact, many compact HII regions as well 
as maser sources have been found from the cloud. 



These sources trace thc^ star forming site in the 
cloud, and are found from the ridge running north 
to south. Since young stellar objects (YSO), re- 
gardless low mass or high mass stars, are rather 
strong X-ray sources with luminosity ranging from 
1032-33 (massive YSO) to lO^s-^s erg s"! (low 
mass YSO) (e.g., Feigelson & Montmerlc 1999), 
numerous unresolved YSOs in Sgr B2, if any, may 
mimic the diffuse X-rays. The spatial distribution 
of the diffuse X-rays, a concave-shape, however is 
largely different from that of the north-south ridge 
of star formation. The X-ray spectra of YSOs are 
generally consistent with optically thin thermal 
plasma of temperature ranging from 1 keV to a few 
keV, with the 6.7 keV line of He-like irons present 
in some cases. The spectrum of the Sgr B2 diffuse 
is completely different from that of typical YSOs, 
hence a YSO origin is ruled out for the bulk of the 
diffuse X-rays. Partial contribution of YSOs may 
still be conceivable in the two local enhancement 
at the cloud center and north, Sgr B2 Main and 
North (see Figure 1). We hence examine the X-ray 
spectrum from these two regions and find no es- 
sential difference from the total diffuse spectrum. 
We only see a hint of a weak 6.7 keV line with the 
equivalent width of 0.10±0.07 keV. Consequently, 
we can safely conclude that no significant contri- 
bution of YSOs is found in the spectrum and mor- 
phology of the diffuse X-rays of Sgr B2. 

4.2. XRN Due to Internal Point Sources 

It may be conceivable that the resolved point 
sources in Sgr B2 irradiate the surrounding gas 
and produce the XRN-like X-rays. These point 
sources, however, are uniformly distributed in the 
cloud (see Figure 1), which disagrees with the 
concave-shape of the diffuse X-rays. The inte- 
grated luminosity of all the resolved point sources 
is 3x 10'^^ erg s^^, which is only 3% of the X-ray 
luminosity of the XRN. Therefore internal point 
sources can be ignored in the XRN scenario. 

4.3. XRN Simulation by an External X- 
ray Source 

The diffuse X-rays are located on the south- 
west half of Sgr B2, the Galactic center side with 
the concave-shape pointing to the Galactic cen- 
ter. These naturally lead us to suspect that an 
irradiating source should be located in the Galac- 
tic center direction. More quantitatively, we nu- 
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merically simulate the XRN model with the as- 
sumption that the gas distribution of Sgr B2 is an 
elliptical cylinder-shape with the same radial de- 
pendence given in the model by Murakami et al. 
(2000) and is irradiated by an external source at 
the Galactic center direction. 

In the XRN scenario, X-rays above the 7.11 keV 
energy are absorbed mainly by iron atoms, while 
the lower energy X-rays are absorbed mainly by 
lighter elements such as carbon, nitrogen and oxy- 
gen. Continuum X-ray flux by the Thomson scat- 
tering is proportional to the number of electrons 
(^hydrogen atoms). Accordingly, the ratio of ab- 
sorption depth at 7.11 keV and low energy band 
is proportional to the abundance ratio of iron and 
the other elements. The equivalent width of the 
6.40 keV line constrains the iron abundance. We 
thus simulate the spectrum with various abun- 
dance sets of iron and the lighter elements and 
compare with the observed spectrum. The best- 
fit spectrum is obtained with the reduced of 
1.14 with 21 degrees of freedom (d.o.f.). The re- 
sult is given in Figure 2 by solid histogram to- 
gether with the data (crosses), where the dashed 
line is contribution of background point sources. 
The abundances of iron and lighter elements are 
2.ltQ g and 3.5^0.? solar, respectively, consistent 
with, but more accurate than the ASCA results 
(Murakami et al. 2000). 

Figure 3 is a simulated XRN image of the 
6.4 keV line with the best-fit abundances. We 
see a clear concave-shape emission, similar to the 
observed result. For quantitative comparison, we 
clearly need more complicated assumption of gas 
distribution, which is beyond the scope of this pa- 
per. 

Using the simulation, we estimate the lumi- 
nosity of an irradiating source should be ^ 10^'' 
(i'/lOOpc)^ erg s~^, in order to fully account for 
the flux of the diffuse X-rays, where £> is a distance 
from the cloud center to the irradiating source. No 
adequately bright X-ray source, however has been 
found near the Sgr B2 cloud (this work) nor at the 
central 20 pc region from Sgr A* (e.g., Baganoff 
et al. 2001). This strengthens the idea proposed 
by Koyama et al. (1996) and Murakami et al. 
(2000) that a past activity of the Galactic center 
is responsible for the presence of reflected X-rays 
and absence of an irradiating source. The past 
activity of Sgr A* would be generated by an ac- 



cretion surge on the MBH due, for example, to the 
passage of dense shell of a young supernova rem- 
nant Sgr A East as is discussed by Maeda et al. 
(2001) based on the new Chandra results of the 
GC. 

5. Summary 

1. Diffuse X-rays are clearly detected from the 
giant molecular cloud Sgr B2, by resolving one and 
half dozens of point sources. 

2. The morphology of the diffuse X-rays is 
concave-shape at the Galactic center side and 
pointing to this direction. 

3. The spectrum exhibits strong emission lines 
at 6.40 and 7.06 keV, a large absorption at low en- 
ergies below 4 keV, equivalent to 9 x 10^^ H cm~^ 
(with the solar abundances) and a pronounced 
edge-structure at 7.11 keV with equivalent TVpe of 
3 xlO^^ cm-2. 

4. The resolved point sources are likely located 
near or in the cloud. The integrated luminosity of 
all the point sources is 3 x 10'^^ erg s~^, which is 
only 3 % of the diffuse X-rays. 

5. The spectrum and morphology of the diffuse 
X-rays can not be attributable to young stellar 
objects in the cloud, nor the fluorescent by the 
resolved point sources. 

6. The Chandra results strengthen the X-ray 
Reflection Nebula (XRN) model, in which a strong 
X-ray source at the Galactic center side irradiates 
the Sgr B2 cloud. 

7. No luminous X-ray source to produce the 
XRN is found near the Sgr B2 cloud. Possible 
idea, as already proposed by Koyama et al. (1996) 
and Murakami et al. (2000), is that the Galactic 
center Sgr A* was active in the recent past but is 
quiescent at present. 
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Fig. 1. Adaptivcly smoothed ACIS-I imago in the 6.15 6.55 kcV band. The Sgr B2 molecular cloud is 
indicated by the contours of the ^'^CO line flux (Sato et al. 2000). The diffuse X-rays are seen in the south- 
west side of Sgr B2 (right side in the figure) with a concave-shape pointing to the Galactic center direction. 
The spectra of the diffuse X-rays and background arc taken from the solid rectangle (3' x 3.'5) and the dotted 
circle, respectively. The positions of the resolved point sources and the data area for the combined spectrum 
are shown by the small solid circles with the radius of six times of the HPD of PSF. 




Energy (keV) 



Fig. 2. — The spectrum of the diffuse X-rays from the Sgr B2 cloud and the best-fit XRN model convolved 
with the response function (solid line). Residuals are shown in the bottom panel. The dashed line indicates 
the contribution of all the resolved point sources in Sgr B2 (see Figure 1). 
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Fig. 3. — The simulated brightness distribution of the 6.4-keV hue overlaid on the eontours of the molecular 
gas model. An irradiating source is assumed to be in the Galactic center direction (the right hand side; 
shown by the arrow). The solid rectangle is the same region as that of Figure 1. 
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Table 1 

Best-Fit Results of the Combined Spectrum of Point Sources 



(10^3 H cm-2) 


(keV) 


Fx' 

(10-5ph cm-2) 


(10^3 erg s-i) 


-■-•^-0.4 


°-0-3.7 


1 c+O.S 
■■■•^-0.2 


^•■'-0.8 


Reduced x'' (d.o.f.) 






1.02 (12) 



Note. — Errors are at 90% confidence level. 
^Hydrogen column density. 
'"Temperature of a thin thermal plasma. 
'^Flux (no correction of absorption) in the 2-10 kcV band. 
Absorption corrected luminosity in the 2-10 keV band 
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Table 2 

Best-Fit Results of the Diffuse X-Ray Spectrum 



Continuum 







pc 


Fx^ 






(10^3 H cm-2) 


(10^^ Fe cm-2) 




(IQ-^ph s"^ cm" 






8.8 +li 


3 4 +3.6 

-2.2 


2.0 (fixed) 


1 9 +0.2 








Iron Lines 






Total 




Ka 














F8 


Fx"^ 






(keV) 


(IQ-^ph s"^ 


(keV) 


(IQ-^ph s"^ cm" 




(10^^ erg s"^) 


qo+0.02 


re (5 +0.5 
-0.6 


'•"^-0.01 


1 1 +°-^ 

-■■•^ -0.1 




1 0+"-'' 

-^•'-'-0.2 


Reduced (d.o.f.) 










1.08 (19) 



Note. Errors are at 90% confidence level. 
"Hydrogen column density. 
''Iron column density. 

'^Photon index of a power-law model for the continuum spectrum. 
''Flux (no correction of absorption) in the 4-10 keV band. 
''Center energy of iron Ka line. 
^Absorption corrected fiux of iron Ka line. 

^Center energy of iron K/3 line. This value is assumed to be 1.103 xFxa (see text). 
''Absorption corrected fiux of iron K/3 line. This value is assumed to be 0.1348 xFkq (see text). 
'Absorption corrected luminosity in the 4-10 keV band. 
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